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Abstract. We study the dynamical appearance of scaling solutions in rela- 
tivistic hydrodynamics. The phase transition effects are included through the 
temperature dependent sound velocity. If a pre-equilibrium transverse flow is 
included in the initial conditions, then it may reach the form of the asymptotic 
Hubble flow, r/t, in short evolution times, 7-15 fm. The numerical solutions are 
found to support to the freeze-out models (Blast- Wave, Buda-Lund, Cracow). 
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1. Introduction 

This article is a contribution to the Quark Matter 2005 poster proceedings, prepared 
with W. Florkowski and T. Csorgo, and based on our recent paper 

The successful parameterizations [0 EH ED of the freeze-out process at RHIC 
indicate that the system formed in Au+Au collisions is highly thermalized and 
exhibits strong longitudinal and transverse expansion that may be characterized by 
the Hubble law, v = Hr. Such form of expansion is used in cosmology, with H 
being a function of time (in spite of the time dependence, H is commonly called the 
Hubble constant). The analogy between the Hubble law used in cosmology and in 
relativistic hydrodynamics is related to the dependence of the four-velocity u M on 
the space-time position of the fluid element 



Here r = ^Jt 2 — x 2 — y 2 — z 2 is the invariant time. The form of the four-velocity 
displayed by Eq. is known as the scaling solution. In many cases it represents 
the asymptotic solution of the hydrodynamic equations [ |3J H3 EI • I n our approach 
we assume that the elementary parton-parton collisions that lead to thermalization 
of the system also lead to collective behavior, i.e., to the development of the pre- 
equilibrium transverse flow. This type of the flow has to be included in the initial 
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conditions for the hydrodynamic equations. This is achieved by the use of the 
following transverse-flow profile at the initial time to = lfm, 

v° r = v r (x,t = to)= H ° r ■ (2) 
y/1 + H 2 r 2 

The quantity Hq may be interpreted here as the Hubble constant that controls the 
amount of the pre-equilibrium flow. In the general case, Hq differs from H which 
is determined by the flow profile at much later times. 



2. Hydrodynamic equations 

Here we briefly introduce the hydrodynamic equations in the way they were formu- 
lated in Ref. The conservation of energy and momentum, d^T^ = 0, in the 
case of an ideal fluid with zero net baryon density has the following form 

u*d tl {Tu v ) = d v T, 9„(aO = 0, (3) 

where T is the temperature, a is the entropy density and — 7(l,v) is the 
hydrodynamic four- velocity (with 7 = l/yl — v 2 ). The role of equation of state is 
played by the temperature dependent sound velocity. 

c 2 (T) = = °™ (4) 

For boost-invariant systems with cylindrical symmetry, Eqs. Q may be put into 
a compact form. In this case: v r = tanha, where a is the transverse rapidity 
and the longitudinal velocity is equal to v z = z/t. With the help of auxiliary 
functions a± = exp ($ ± a) and the $ potential defined by the differential formula 
rf$ = d In T/c s = c s d\n<r 1 the characteristic form of the hydrodynamic equations is 
introduced 

|a ± (t, r) + 1 a ± (t, r) + ( ^ + I) a ± (t, r) = 0. (5) 

at 1 ± v r c s or 1 ± v r c s \ r t J 

Knowing the solutions of this equation, i.e., the functions a±, we are able to calculate 
the transverse velocity and the $ function from the expressions 

$ = iln( a+ a_). (6) 



a + + a_ 2 



2.1. Temperature dependent sound velocity 

We generalize the approach introduced in [ |H] by the use of the temperature de- 
pendent sound velocity in the numerical calculations In order to evaluate the 
temperature from the auxiliary functions a± we first obtain $ as a function of the 
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temperature. We denote this function as $y With the help of the inverse function, 
denoted by T$, we calculate the temperature from the formula 

T (r, i) = T»Q In (a+o)- (7) 

In our numerical calculations we consider two forms of the sound velocity. The first 
one is obtained from the QCD lattice calculations [ |5] and extrapolated to lower 
temperatures, while the second one is the analytic parameterization of these points. 
The latter form allowed us to calculate functions <&t and Xb analytically. 

2.2. Initial conditions 

We use the same method for constructing initial conditions as in i.e., a single 
function a(r) is used to determine the initial forms of the functions a+ and a_, 
namely a+(t = to,r) = a(r), a_ (t — to,r) — a(— r). Such procedure insures that a 
cylindrically symmetrical system has vanishing velocity field at r = 0. In addition, 
we assume that the initial temperature profile is connected with the thickness func- 
tion Ta (r) . We assume that the initially produced entropy density is proportional 
to the number of nucleons participating in the collision at the distance r from the 
collision center, cr(r) ~ Ta(t), and at high temperature our system may be consid- 
ered as a system of massless particles. Since the entropy density of a massless gas 
is proportional to T 3 we find 

T(r,t = t ) = const Tf (r) . (8) 
The final form of the function a(r) that satisfies initial conditions @ and JSJ reads 



(9) 

3. Results 



a(r,t = t Q ) = a T (r) 



V 1 
. PC 



ar(r) = cxp 



$ T (const Tf (|r|) 



Our numerical studies show that the absence of the initial pre-equillibrium flow, 
the case Ho — 0, unables formation of the scaling solutions during short evolution 
times, t < 10 fm (see the upper parts of Fig. [IJ. In the case where the pre- 
equilibrium flow is included in the initial conditions we observe that the velocity 
profiles may approach the asymptotic scaling solutions after about 7 fm (see the 
lower parts of Fig. The temperature profiles obtained in all considered cases 
show similar behavior; the system cools fast until it reaches the critical temperature 
and further cooling is slowed down. However, in the case of the non-zero initial flow 
the decrease of the temperature is faster and after some time the inner part of the 
system has lower temperature then the edge. This leads to the hyperbolic shape 
of the isotherms, as shown in the lower part of Fig. ^b). We conclude that this 
type of behavior gives support for the parameterizations used in the Blast- Wave [ 
E], Buda-Lund [0, and Cracow [0] models. 
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Fig. 1. Two cases of the hydrodynamic expansion of matter. The first case corresponds 
to matter being initially at rest, Ho = 0, and the second case includes the pre-equilibrium 
flow characterized by the velocity profile (2) with Ho — 0.25 fm -1 . The initial central 
temperature To equals 2T C (T c — 170 MeV), while the sound velocity is taken from the 
QCD lattice calculations. The parts a) show the transverse velocity as functions of the 
distance from the center for 6 different values of time, ti = 1,4, 7, 10, 13 and 16 fm. The 
dashed thin lines describe the ideal Hubble-like profiles of the form r/U (r < ti). The parts 
b) show the temperature profiles in r, whereas the parts c) show the isotherms in the t — r 
space. In this case, the labels at the curves denote the values of the temperature in units 
of the critical temperature. 
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